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Topological superconductors may undergo transitions between phases with different topological 
numbers which, like the case of topological insulators, are related to the presence of gapless (Ma- 
jorana) edge states. In Z topological insulators the charge Hall conductivity is quantized, being 
proportional to the number of gapless states running at the edge. In a superconductor, however, 
charge is not conserved and, therefore, a^y is not quantized, even in the case of a Z topological 
superconductor. Here it is shown that while the axy evolves continuously between different topo- 
logical phases of a Z topological superconductor, its derivatives display sharp features signaling 
the topological transitions. We consider in detail the case of a triplet superconductor with p-wave 
symmetry in the presence of Rashba spin-orbit (SO) coupling and externally applied Zeeman spin 
splitting. Generalization to the cases where the pairing vector is not aligned with that of the SO 
coupling is given. We generalize also to the cases where the normal system is already topologically 
non-trivial. 



PACS numbers: 71.70.Ej, 74.25.fc, 73.43.-f 
I. INTRODUCTION 

In two dimensions, Z topological insulators exhibit a 
charge Hall conductivity that is quantized and propor- 
tional to the Chern number of the occupied bandai"^. 
Such nontrivial topological phases are a lso characterized 
by the presence of gapless edge mode^^ that can be 
detected by transport measurements or tunneling. In a 
topologically non-trivial superconductor one does not ex- 
pect that the charge Hall conductivity may be quantized, 
however, as charge is not conserved due to the breaking 
of C/(l) symmetry. In a singlet superconductor spin is 
conserved and there is still a possibility that the spin 
Hall conductivity is quantized, as previously shown for a 
d-wave superconductor in the vortex state^. For a triplet 
superconductor even this quantization is absent. The 
thermal Hall conductivity has recently been shown to be 
quantized for topological superconductors with broken 
time reversal symmetry (TRS^. 

Generically speaking, the charge Hall resistance may 
be written as the sum of two contributions, one pro- 
portional to the magnetic field and an anomalous con- 
tribution as pxy = RqHz -\- (considering z as the 
perpendicular direction to the plane where the charges 
move) . The term p^y is the anomalous Hall effeclP and 
has different origins. One of these origins is intrinsic 
and is the result of an anomalous velocitjl^ that is the 
result of a non-zero Berry curvatur^l^, r2„(fc). The ve- 
locity of a charged particle in a given energy band n 
in the presence of an electric field, E, can be written 
as t;„(fc) = h-^de„{k)/dk - {e/h)E x rj„(fc). The last 
term gives a contribution to the velocity that is trans- 
verse to the direction of the electric field and therefore 
contributes to the Hall conductivity. Two other mecha- 
nisms that lead to an anomalous velocity are due to scat- 
tering from impurities in a system where the spin-orbit 
has to be taken into account such as the skew scattering 



mechanisrrPi' and the side jumfP^. 

Historically, the anomalous Hall resistivity was stud- 
ied in detail in systems with a finite magnetization Af^, 
where p^^ = i?sM2.'^ Anomalous properties in supercon- 
ductors with magnetization have been studied before in 
particular, the presence of magnetoelectric eSect^^, the 
generation of a charge Hall effect due to a spin currenlP^ 
or a spin Hall effect due to a charge currentP^. The 
anomalous Hall effect has been studiecP^ in supercon- 
ductors with spin-orbit coupling. The presence of a mag- 
netic impurity is enough to induce a non-vanishing Hall 
conductivity^^. The magnetic impurity interacts with 
the superconductor as a local Zeeman term and orbital 
effects such as the presence of vortices are neglected. The 
local magnetization may also be the result of some prox- 
imity effect for instance with a magnetic dot. A dense 
magnetic impurity distribution leads to the destruction 
of superconductivity if the pairing is a spin singlet but 
magnetization and superconductivity may coexist if the 
pairing is of triplet origin such as in a p-wave supercon- 
ductor. In this case the Hall conductivity was shown to 
be non-vanishing if the spin-orbit coupling is presenlPSl. 

In this work we will be concerned with the effect of 
the intrinsic contribution to the charge Hall conductiv- 
ity in Z topological superconductors. Superconductiv- 
ity with non-trivial topology may be obtained in differ- 
ent ways.l^ It can be due to the pairing symmetry, as is 
the case of p-wave superconductors.'^ In semiconductors 
with Rashba spin-orbit coupling it arises when s-wave su- 
perconductivity is induc ed an d a Zeeman (time-reversal 
breaking) term is added.'^*'^^' An interesting proposal is 
that of systems where the normal phase is already topo- 
logically non-trivial, in which case a topological super- 
conductor can be obtained if s-wav e superconductivity 
is induced by proximity effect .'^^'^ Here we consider a 
Rashba-type non-centrosymmetric superconductor with 
admixture of s-wave and p-wave pairing and TRS break- 
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ing Zeeman splitting, which has been recently proposed 
in Ref. [HI We generalize to the cases where the pairing 
vector is not coincident with that of the spin-orbit cou- 
pling, and also to the cases where the normal phase is 
topologically non-trivial. We thus take on equal footing 
the three possible ways of obtaining topological super- 
conductivity mentioned above. 

The main results of this paper may be summarized as 
follows. Whenever a Z topological superconductor is re- 
alized and the first Chern number fully characterizes the 
topological phase, we find that the behavior of the Hall 
conductivity and its derivatives with respect to the pa- 
rameters that drive the topological phase transition, spe- 
cially the second derivative, provide an alternative way to 
identify topological transitions in superconductors. This 
approach proves extremely useful when the pairing vector 
is not aligned with the spin-orbit, in which case we have 
found a less obvious relation between the Chern number 
and the number of crossings of edge state bands with the 
Fermi level. A careful topological analysis of this case is 
given. 

The paper is organized as follows. In Section II we 
introduce the expression for the Chern number that will 
be used throughout this work. The model Hamiltonian 
is presented in Section III and in Section IV the re- 
sults for the Hall conductivity and Chern number are 
presented. Section V is devoted to the analysis of the 
edge states and their correspondence to the topological 
indices. In Section VI we present results for a model 
where the non-superconducting band structure is already 
nontrivial. Our conclusions are presented in Section VII. 



II. CHARACTERIZATION OF TOPOLOGICAL 
PHASES 

The Berry curvature tensor for a band with Bloch 
wavefunctions Un{k) can be calculated as 



J^n(fc) = (VfeM„(fc)| X |V 



(1) 



where k = {kx,ky) denotes the momentum vector. The 
contribution from the n-th band to the Hall conductivity 
of a normal system may be written in terms of the Berry 
curvature aa^ 
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where np is the Fermi function. If the chemical poten- 
tial lies within a gap the integral over the occupied states 
runs over the entire Brillouin zone. The charge Hall con- 
ductivity can then be written as 



(3) 



where C is the sum of the Chern numbers of the occupied 
bands. The Berry curvature may also be obtained as a 



sum over states analogous to the Kubo formula for the 
conductivity, and reads: 
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In the case of a superconductor, the states \n) are the 
eigenstates of the Bogoliubov-de Gennes equations. At 
the gapless points the denominator vanishes and the in- 
tegral over the Brillouin zone may have large numerical 
errors. It is then more convenient to calculate the Chern 
number by computing the flux of the Berry curvature 
over plaquetes in the Brillouin zon^l^. Discretizing the 
Brillouin zone as = iirj/N, with j = l,...,iV, and 
fi ~ x,y, a new variable, U^{k)^ for the link 5k^ (with 
l^fc^l = 2tt/N) oriented along the ji direction from the 
point k may be defined as 



(n(fc)|n(fc + Jfc^)) 



\{n{k)\n{k + 5k^))\ 
and the lattice field strength may be defined as 



(5) 



F^y{k) = In {Ux{k)Uy{k + 5k^)Ux{k + 5ky)-^Uy{k)-^) . 

(6) 

Fr^yik) is restricted to the interval — tt < ~iFxy{k) < tt 
and the gauge invariant expression for the Chern number 
is 



(7) 



The calculations of the Chern number of each band n are 
performed in this way in this work. 

The calculation of the Chern number is simple for a 
2x2 Hamiltonian matrix H once the latter is written in 
the form: 



H{h) = h{k) 



ho{k)To , 



(8) 



where h = {hx, hy,hz), t are Pauli matrices and tq is the 
identity. The Chern number for the bands in Hamilto- 
nian Eq. Q is independent of the choice for ho{k), as 
computed from the usual expression 
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dh ^ dh 



(9) 



h — h/\h\. The topological nature of bands may be 
understood as the result of the covering of the unit sphere 
defined by the vector h. 

The system's symmetry properties depend on whether 
the Pauli matrices in equation ([s]) represent a pseudospin 
{e.g., a sublattice) or the physical spin. In the first case 
TRS requires hx(^z) to be an even function of k and hy 
to be odd. Otherwise, all components have to be odd. 
In order to have nonzero C, TRS must be broken. The 
operation of spatial inversion does not change the Chern 
number C. 
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On general grounds, non-trivial topological order for 
non-interacting Hamiltonians can be related with the 
presence or absence of three discrete symmetries: time- 
reversal, particle-hole, and chiral symmetryP^l^^l For 
Bogoliubov-de Gennes systems, where particle-hole sym- 
metry is always present, preserving or not TRS is de- 
terminant to the nature of possible topological phases 
in two dimensions. The non-centrosymmetric supercon- 
ductor we consider here is time-reversal invariant if the 
Zeeman term is absent and the pairing is unitary. The 
system then belongs to the symmetry class Dili where 
the topological invariant is a Z2 index, and it is said to 
realize a Z2 topological superconductor. If the pairing is 
non-unitary or the Zeeman term is finite TRS is broken 
and the system belongs to the symmetry class D (the 
TRS operator T is such that = —1). The topological 
invariant that characterizes this phase is the first Chern 
number C, and the system is said to be a Z topological 
superconductor. 



III. TOPOLOGICAL SUPERCONDUCTOR 

We consider a triplet superconductor with p-wave sym- 
metry in the presence of Rashba spin-orbit coupling and 
magnetization, e.g., due to a time- reversal breaking Zee- 
man term. Due to the non-centrosymmetric nature of 
the system, parity is broken and, in general, the pair- 
ing symmetry is not fixed, and an admixture of singlet 
pairing is allowed^^. Therefore, we also consider a con- 
tributi on fro m s-wave pairing. This model was studied 
in Refs.'i^'^S] -^j-jte the Hamiltonian as 



H 



-E 

k 
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Here, ek — —2t{coskx + cosky) — ep is the kinetic part, 
where t denotes the hopping parameter set in the follow- 
ing as the energy scale, t = 1, cp is the chemical poten- 
tial, fc is a wave vector in the xy plane, and we have taken 
the lattice constant to be unity, a = 1. Furthermore, 
is the Zeeman splitting term responsible for the magne- 
tization, in energy units, along the z direction. Finally, 
the Rashba spin-orbit term is written as 

Hp; = s ■ (J = a (sinfcyCTj; — smk^Oy) , (12) 

where a is measured in the energy units and s ~ 
Q!(sin fcj^, — sin fcj,, 0). The matrices ctj; , CTj, , cr^ are the 
Pauli matrices acting on the spin sector, and uo is the 
2x2 identity. 

The pairing matrix reads 



A = « (d • cr + As) (Jy 



rf. - A 



'A, 
id,, 



(13) 



The vector d = (dx, dy,dz) is the vector representation of 
the p-wave superconducting pairing and is an odd func- 
tion of fc. Because of Fermi statistics, the pairing matrix 
satisfies A(fc) — — A-^(— fc). The triplet pairing term is 
invariant under a spin rotation about the d direction. We 
note that both the superconducting order parameter and 
the magnetization may be due to intrinsic order or to 
some proximity effect due to neighboring superconduc- 
tors or ferromagnets. The pairing matrix for a p-wave 
superconductor generally satisfies 



AA^ = |d|2 



(To + q • T , 



(14) 



where q — id x d* . If the vector q vanishes the pairing 
is called unitary (s-wave pairing is always unitary) . Oth- 
erwise it is called non-unitary2Sl and breaks TRS, orig- 
inating a spontaneous magnetization in the system due 
to the symmetry of the pairing, as in '^He. 

We will consider both unitary and non-unitary pair- 
ings. In the case of unitary pairing we consider two ex- 
amples. One of them respects to a situation where the 
spin-orbit coupling is such that the pairing is alignecP^ 
along the spin-orbit vector s. This is a situation expected 
if the spin-orbit is strong since it is energetically favor- 
able, and we will denote it by strong coupling case. In 
the other case we will relax this restriction and allow that 
the two vectors are not aligned. This case we will denote 
by weak spin-orbit coupling. 

The energy eigenvalues and eigenfunction may be ob- 
tained solving the Bogoliubov-de Gennes equations 
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The 4-component spinor can be written ai 



mi 
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w«(fc,t) 

Un(k, \,) 

-w„(-fc,t) 

Vn{-k,i) 



(16) 



The energy eigenvalues of the Hamiltonian ( 10 ) can be 
written (for As = and d^ — 0) as 



where 

zi 



i-k,i 



= dd* 



(17) 



s- s + ei+ Mt 



= \{dxs)z\ - iekMz{dx d*)z 



[d ■ d*)2 - \d ■ dr + e|(s • s + Ml), (18) 



and Qfi, a2 = ±- The gap between the lowest bands closes 
at the fc points satisfying the condition zi = 2,Jz2- 

In the superconducting phase the system is generally 
gapped. A possible change of topology occurs when the 
gap closes. Considering the strong spin orbit case for 
which the d and s are collineai^, we may write d = 
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FIG. 1: (color online). Hall conductivity in a unitary case 
with dx — d sin ky,dy = d sin kx,dz — (top panel) and a non- 
unitary case with dx = — d/2sinfca; = —idy,dz = (bottom 
panel) with d= l,eF = —1- This Figure corrects a previously 
obtained resullP^. 

{d/a)s. Taking also the s-wave pairing into account, the 
gapless points satisfy 

€k-s = A,s. (19) 
a 



IV. CHERN NUMBERS AND HALL 
CONDUCTIVITY 

A. Weak spin-orbit coupling: Unitary and 
non-unitary pairings 

In Fig. [l] we present the results for the Hall con- 
ductivity, as calculated from the Kubo formula, given 
in Eq. (10) of Ref.^^. We consider both the unitary 
and the non-unitary cases, relaxing the restriction that 
d II s which was assumed in Ref. [18l The two cases 
are chosen as = d sin ky,dy — d sin k^,dz = and 
dx = —d/2smkx = —idy,dz = 0, respectively. The Hall 
conductivity in the superconducting unitary phase is sim- 
ilar to that of the normal phase (not shown). Both in the 
normal phase and in the unitary case a magnetization is 
required in order to have a non- vanishing Hall conductiv- 
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FIG. 2: Chern number for the unitary case of Fig. [T] as a 
function of magnetization for d = 1, ep = —1 and a = 1. 

ity. However, in the non-unitary case, the magnetization 
induced by the vector q produces a finite Hall conductiv- 
ity even if the explicit Zeeman term is absent. In all cases 
the spin-orbit coupling is necessary for a non-vanishing 
Hall conductivity. 

In all three cases (normal, unitary, and non-unitary) 
the Hall conductivity has a clear minimum when the 
magnetization is of the order of the chemical potential. 
At this point the spectrum is gapless and, as shown in 
Fig. [2] for the unitary case, the Chern number of the 
occupied bands changes. The topological transition does 
not depend on a and the results in Fig. [2] are shown 
for the particular case of a = 1. It turns out that 
for the non-unitary case depicted in Fig. [T] the spec- 
trum is always gapless due to the lack of dependence 
of the pairing function on ky. Therefore, the expres- 
sion for the Chern number suffers from numerical in- 
stability. Considering a non-unitary pairing of the form 
dx = d sin ky,dy ~ id sin kx, dz = the Hall conductivity 
is similar to that obtained in Fig. [l] and the change in 
the Chern number is also similar to the one shown in Fig. 
[2] for the unitary case. 

In the normal phase the system is topologically trivial 
and the Chern number is zero throughout the space of 
parameters of the chemical potential and the magnetiza- 
tion. 



B. Strong spin-orbit coupling 

For strong spin-orbit coupling it has been shown that 
it is more favorable that the pairing vector d aligns with 
the spin-orbit coupling: 

d{k) ^ d{sin ky, — sin kx) ■ (20) 

As a consequence, the critical temperature associ- 
ated with this type of pairing is higheif^. There is 
a rich sequence of topological transitions as a func- 
tion of the chemical potential, spin-orbit coupling and 




M =0.5 



magnetizatiorflSl. In general this problem involves solv- 
ing for the eigenvalues of the 4x4 matrix in Eq. ( 10 1 



The calculation of the Chern number of each band is 
performed using the eigenfunctions of this 4x4 matrix 
m Eq. ([Tf. Since the gap must close at the topologi- 
cal transitions, the location of these transitions may be 
determined looking at the gapless k point&i^ satisfying 
Eq. (19 1. The location of the transitions and the associ- 



ated gapless points in the spectrum have been obtained 
before^®. It turns out that in each topological phase, the 
Hamiltonian can be continuously deformed in such a way 
that (a, As — ?> 0) without closing the gap. The problem 
then simplifies since the original 4x4 Hamiltonian has 
been deformed to two 2x2 matrices and the Chern num- 
ber may be calculated as in Eq. ^. 

In Fig. [3] we show the results for the Chern number 
of the occupied bands as a function of the chemical po- 
tential and magnetization. There are various transitions 
that correspond to the closing and opening of gaps in 
the spectrum. In the same figure we show the results 
for the Hall conductivity for the same region of parame- 
ters. Even though the Hall conductivity is a continuous, 
smooth function there are clearly local maxima and min- 
ima that can be associated to the points where a topo- 
logical transition occurs. 

In the case of a Z topological insulator, a topological 
transition modifies the Chern number and the value of 




FIG. 4: (color online). Chern number. Hall conductivity and 
its first and second derivatives as a function of magnetization 
and as a function of chemical potential for strong spin-orbit 
coupling. The Hall conductivity and its first derivative are 
multiplied by a factor of 10 for better visualization. 



the Hall conductivity which, therefore, exhibits a clear 
signature of the transition. In the case of a Z topolog- 
ical superconductor there is, however, no discontinuity 
of a^y but its second derivative signals the transitions 
sharply. This is well illustrated in Fig. [4] where we show 
cuts at constant chemical potential as a function of the 
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FIG. 5: (color online). Chern number and Hall conductiv- 
ity for s-wave pairing and zero p-wave pairing (d = 0, As = 
0.5, a = 0.6). 



case of a topological insulator, the derivative is a Dirac 
delta function and the second derivative is the derivative 
of a Dirac delta function. In the superconductor these 
delta functions are smeared but are still clear evidence 
of the location of the transition and, moreover, of the 
change in the topological number. 

We may also consider the case where there is only s- 
wave pairing and no p-wave pairing. The results for the 
Chern number and the Hall conductivity for the same 
region of parameters are shown in Fig. [Sj As shown be- 
fore, if the magnetization vanishes, and the s-wave com- 
ponent is larger than the p-wave component, the phase is 
topologically triviahfi^ There are non-trivial phases that 
arise due to the presence of the magnetization and the 
consequent breaking of TRS.'^-'^^' Note that there is a fi- 
nite region around zero magnetization in which the Chern 
number vanishes, as mentioned. As above, the Hall con- 
ductivity clearly shows the location of the various transi- 
tions between the Chern numbers and an analysis similar 
to the one carried out for the p-wave case both for the 
Hall conductivity and its derivatives may be carried out, 
signaling in a similar way the various transitions. The 
same holds for weak spin-orbit coupling. The topological 
transitions are also clearly detected by the derivatives of 
the Hall conductivity. 



EDGE STATES 



magnetization or fixing the magnetization and changing 
the chemical potential. The results for the Chern num- 
ber clearly indicate the topological transitions either as 
a function of the magnetization or chemical potential. 
The behavior of the Hall conductivity correlates with 
these transitions. As expected from general considera- 
tions if a transition occurs between Chern numbers of 
different signs, the Hall conductivity changes sign accord- 
ingly. Here we are interested in finding a signature of the 
change in the Chern number. This can be achieved by 
looking at the derivatives of the Hall conductivity. At 
the transitions the derivative behaves in a way qualita- 
tively similar to the case of a Z topological insulator: if 
the Chern number increases across the transition the first 
derivative of the Hall conductivity is positive and if the 
Chern number decreases the derivative is negative. The 
change is small and the features in the first derivative are 
also small, therefore, we have multiplied the Hall conduc- 
tivity and its first derivative by a factor of 10. A much 
stronger signal is provided by the second derivative. At a 
transition where the Chern number changes, the second 
derivative exhibits two close peaks: a negative peak fol- 
lowed by a positive one when the Chern number increases 
(and vice- versa if the Chern number decreases). In the 



Due to the bulk-edge correspondence, complementary 
information on the topological phases and transitions 
may be obtained by analyzing the edge states. We con- 
sider a strip geometry of transversal width Ny and apply 
periodic boundary conditions along the longitudinal di- 
rection, X. We write 




(21) 



and rewrite the Hamiltonian matrix in terms of the op- 
erators (21 1 as 



jy 



t 



The operator H^^^j^ reads 



(22) 
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eF~t7]+ ia sin kx + §ir]^ 

" rj- —2tcoskx + Mz — e.p ^ 

i_ -A, 

idsinfcj, + -§1^1- 



—idsinkx 
tr]+ -As 

2t cos kx + Mz + ep 
ia sin kx + ^rj. 



r 



A 

—id sin ka 
—ia sin fc^ 

2t cos k. 



M 



where V']^77±'/'j„ = V'j^V'j^+i ± The diago- 

nahzation of this Hamiltonian involves the solution of a 
^Ny X A.'Ny eigenvalue problem. The energy states include 
states in the bulk and states along the edges. 



A. Strong spin-orbit coupling 

In the case of strong spin-orbit coupling with — 
there is no TRS bre aking and the system belongs to the 
symmetry class 0111.123221 Jn the s-wave case there is only 
the bulk gap and no gapless (edge) states. The system 
is in a topologically trivial phase. In the case of p-wave 
pairing even though the Chern number vanishes there are 
gapless edge stateJ^. The system is in a Z2 topological 
phase. The gapless edge states have a twofold Kramers 
degeneracy and two counterpropagating edge modes give 
opposite contributions to the total Chern number, C = 0. 
This is a similar situation to that in the spin Hall effect, 
where, even though the charge current vanishes, there is 
a spin current along the edges. In the case where there 
is a mixture of s- and p-wave components and the ampli- 
tude of the p-wave pairing is larger than the correspond- 
ing amplitude of the s-wave case, there are edge states 
and a topologically nontrivial phase. Because of spin- 
momentum locking there is no backscattering and these 
states are topologically protected from non-magnetic im- 
purities. 

As the magnetization is turned on TR S is br oken and 
the system's symmetry class changes to Dl ^^l^^l For small 
magnetization the Z topological superconductor is in a 
trivial phase with Chern number C = 0. A finite mag- 
netization is then necessary to cause a topologicalphase 
transition to a phase with non-zero Chern number.U^This 
happens both for the p-wave case and the s-wave case. 
The sequence of Chern numbers is clearly correlated with 
the number of pairs of edge states as shown in Ref. UHl 

In Fig. [6]we show the edge modes as a function of spin- 
orbit coupling for the case of strong spin orbit coupling 
(d II s) and small magnetization, when the system is in 
the C = phase. The spin orbit coupling does not change 
the Chern number since it does not close the bulk gap. It 
is interesting to note that even though the system is in a 
C = phase the number of edge states is two; the same as 
that in the parent Z2 phase, when = 0.1^ For 7^ 0, 
however, TRS is broken and these edge states are not 
topologically protected against (any type of) disorder. In 
this sense the system is in a trivial phase, in accordance 
with the Chern number C = 0. Nevertheless, in the clean 
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FIG. 6: Gapless edge modes for unitary case of strong spin- 
orbit coupling for zero Chern number for different values of 
the spin-orbit coupling, a. Here = —3, d — 0.6, A = 
0.1, = 0.5. 



limit these edge modes could be detected. 

The presence of edge modes induced by bulk topol- 
ogy can also be shown using dimensional reduction and 
thereby calculating the winding number.l^ For fcy = or 
TT, the Hamiltonian H{k) has the chiral symmetry: 



ri/(fe)r^ = -H{k) 



(24) 



with r = Tj; (g) (Toj where ctq is the identity in spin space 
and Tx acts on the particle-hole space. The operator that 
diagonalizes F is^SiMi 



T = a-Q^ e 



(25) 



and the Hamiltonian can then be brought to the off- 
diagonal form: 



TH{k)T^ 



q{k) 



(26) 



if ky — 0,7T and d^ — where 



q{k) 



— — Mz + id sin kx ia sin kx — A^ 

—ia sin kx + A^ efc — ef + + id sin k. 



(27) 



8 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 
k/jr k/n k/ji 



0.5 j 

o| 

■0.5 1 

0.6i 
0.31 

bo 

u 0| 

w -0.3I 

-0.61 



0=0.1 




J y 




A h 





-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 
k/JC k/7C k/3t 




FIG. 7: Gapless edge modes in the unitary case for zero Chern 
number (top) and Chern number equal to 2 (bottom). Here 
= -l,d = 1 and = 0.5 (top), = 1.2 (bottom). 



FIG. 8: (color online). Gapless edge modes for unitary case 
for zero Chern number (top) and Chern number equal to 2 
(bottom). Here ep = — l,d = 1 and = 0.5, a = 2 (top), 
and = 1.2, a = 3 (bottom). 



The winding number is then defined as 



I{ky) 



(ky = 0,7r) 



(28) 

Physically, a nonzero I{ky) means that if the system is 
infinite along the y direction and finite along x, there 
will be edge states with ky — or The calculation 
of the winding number gives the number of gapless edge 
modes both when the Chern number vanishes and when 
the Chern number is finite.'i^ 



B. Weak spin-orbit coupling 

If the spin-orbit coupling is not strong, so that the 
pairing vector d is not aligned with the spin-orbit vec- 
tor, as in the unitary and non-unitary cases considered 
in Sec. |IV[ the connection between the Chern number 
and the number of gapless edge states is less transpar- 
ent. Fig. [7] shows the low- lying energy modes for the 
unitary case previously considered for different values of 
the spin orbit coupling, a. The top panel corresponds to 



a case where the Chern number vanishes while the bot- 
tom panel to a non-vanishing Chern number. There is 
a variety of gapless edge states for both the C = and 
C = 2 cases. Even though changing a should not change 
the topology, there is an apparent appearance of vari- 
ous gapless states that seem not to follow the bulk-edge 
correspondence. However, in the unitary case Eq. (25) 



still transforms the Hamiltonian into an off-diagonal form 
similar to that in Eq. (27), so that the winding number 
is still well defined. A calculation of the winding number 
shows that the number of gapless edge modes is actually 
independent of a. In the top panel (C = 0) we get that 
/(O) = 2 and /(tt) = and in the bottom panel we obtain 
that /(O) = 1 and /(tt) = 1 in agreement with the value 
C = 2. 

The bulk-edge correspondence is further elucidated in 
Fig. [8] Careful analysis shows that some of the gapless 
edge states do not originate from bulk topology and that 
the number of topologically induced edge states (given 
either by the winding number or the Chern number in 
the case of C 7^ 0) is consistent. Only the bands of edge 
states that connect the upper and lower bulk bands, i.e. 
connecting open and filled circles in Fig.|8j can be traced 
back to the nontrivial bulk topology.!^ Denoting the two 
edges of the system as R and L we see that for C = 0, 
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FIG. 9: (color online). Chern number and Hall conductivity 
in the case where the normal system is topologically nontriv- 
ial. The parameters are: ti — l,t2 = 1.1, As = 0.1, d = 
0.6, a = 0.6. In the lower panel ep = 1.92. 
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FIG. 10: Gapless edge modes for the same system as in Figure 
|9]with ep = 1.92 and: \h = 1, C = -3 (top); A/^ = 4, 
C = -1 (bottom). 



the number of propagating states at each edge is always 
the same as the number of counter propagating ones. 
For C = 2, on the other hand, the difference between 
propagating and counter propagating states is always 2 
at each edge. 



VI. NONTRIVIAL TOPOLOGY IN NORMAL 
PHASE 

The nontrivial topology of the bands in the supercon- 
ducting phases above stems from the mixture of the par- 
ticle and hole excitations since the normal phase is non- 
topological. We may as well consider a system that is 
nontrivial in the normal phase and add superconductiv- 
ity either self-consistently or via a proximity effect. This 
has been proposed before in various contextJ^^. 

Considering a Hamiltonian of the form of Eq. (|8| and 
selecting 



a sin ky 



-a sin kx , 



hz = 2ti (cos kx + cos ky) + At2 cos kx cos ky , (29) 

leads to nontrivial phases as the hoppings ti and t2 are 
variecpfil. Results for the edge states, the Chern num- 
ber and the Hall conductivity are shown in Figs. [9] 
and 10 In the normal phase TRS is broken since 



is even in the momentum. The Chern number is C = 2 
if < \t2\] and C = 1 if > 1^2 1- As the system be- 
comes superconducting the nontrivial topology remains 
even though the Chern number changes. The nontrivial 
topology of the normal state bands lends some robust- 
ness to the topological superconducting phase. Indeed, 
the Chern number remains invariant in large portions of 
the parameter space. In the first panel of Fig. |9]we show 
cuts of the Chern number at constant chemical poten- 
tial, ep = —1.92,0,1.92, as a function of M^. For neg- 
ative chemical potential the Chern number is C = —3 
except for some narrow regions where C = 5. For zero 
and positive chemical potential there is a single topolog- 
ical transition from C = — 3toC = — 1. The results for 
the edge states in Fig. [TO] show a clear correspondence 
between the Chern number and the number of gapless 
modes. Note that in this case there is strong spin-orbit 
coupling. The difference to the previous Sections is the 
non-trivial topology of the normal phase. The supercon- 
ducting pairing then changes the topology, as shown by 
the change of Chern number entering the superconduct- 
ing phase. In this case the Hall conductivity also signals 
the transition at positive values of the chemical potential 
and varies smoothly in the narrow region where C = 5. 
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VII. CONCLUSIONS 

We have shown that the Hall conductivity and its 
derivatives may be used to detect the topological tran- 
sitions that occur in Z topological superconductors. In 
a topological insulator the Hall conductivity is quantized 
and proportional to the Chern number and, therefore, its 
discontinuous changes across a transition can be used to 
detect and characterize the transition. Even though the 
Hall conductivity is not quantized in a superconductoi^ 
it may also be used to study these transitions. This pro- 
vides a bulk detection method of these transitions that is 



complementary to the detection of the gapless edge states 
associated with these nontrivial topological phases. 

In the case of strong spin-orbit coupling there is a sim- 
ple correspondence between the number of gapless edge 
states and the Chern number, both for trivial and non- 
trivial normal state bands. However, in the case of weak 
spin orbit coupling where the pairing vector, d, is not 
parallel to the spin-orbit vector, s, extra unprotected 
gapless modes appear. The bulk-edge correspondence is 
preserved as evidenced by the calculation of the winding 
number. 
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